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<p ' Abstract 

Prion diseases are invariably fatal and highly infectious neurodegenerative diseases af- 

Qh! footing humans and animals. The neurodegenerative diseases such as Creutzfeldt-Jakob 

disease, variant Creutzfeldt-Jakob diseases, Gerstmann-Straussler-Scheinker syndrome, 
Fatal Familial Insomnia, Kuru in humans, scrapie in sheep, bovine spongiform en- 
cephalopathy (or mad-cow disease) and chronic wasting disease in cattle belong to 



X 



> 

£T} ', prion diseases. By now there have not been some effective therapeutic approaches to 

treat all these prion diseases. In 2008, canine mammals including dogs (canis familials) 
were the first time academically reported to be resistant to prion diseases (Vaccine 
26: 2601-2614 (2008)). Rabbits are the mammalian species known to be resistant to 
infection from prion diseases from other species (Journal of Virology 77: 2003-2009 
(2003)). Horses were reported to be resistant to prion diseases too (Proceedings of 
the National Academy of Sciences USA 107: 19808-19813 (2010)). By now all the 
NMR structures of dog, rabbit and horse prion proteins had been released into protein 
data bank respectively in 2005, 2007 and 2010 (Proceedings of the National Academy 
of Sciences USA 102: 640-645 (2005), Journal of Biomolecular NMR 38:181 (2007), 
Journal of Molecular Biology 400: 121-128 (2010)). Thus, at this moment it is very 
worth studying the NMR molecular structures of horse, dog and rabbit prion proteins 
to obtain insights into their immunity prion diseases. 

The author found that dog and horse prion proteins have stable molecular dynami- 
cal structures whether under neutral or low pH environments, but rabbit prion protein 
has stable molecular dynamical structures only under neutral pH environment. Under 
low pH environment, the stable a-helical molecular structures of rabbit prion protein 
collapse into /3-sheet structures. This article focuses the studies on rabbit prion protein 
(within its C-terminal NMR, Homology and X-ray molecular structured region RaPrP c 
(120-230)), compared with human and mouse prion proteins (HuPrP (125-228) and 
MoPrP c (124-226) respectively). The author finds that some salt bridges contribute 
to the structural stability of rabbit prion protein under neutral pH environment. 



Key Words: Neurodegenerative Diseases; Prion Diseases; Immunity; Rabbit Prion 
Protein, Rabbit NMR, X-ray and Homology Structures. 

1 Introduction 



As we all know, the disease infectious prions PrP bc are rich in /3-sheets (about 43% 



o 



/3-sheet) [Griffith 1967] and the normal cellular prions PrP c are predominant in en- 
helices (42% a-helix, 3% /3-sheet) [Pan et al. 1993| . Prion diseases are believed caused 
by the conversion of normal PrP to abnormally folded PrP , and prion diseases are 
so-called protein 'structural conformational' diseases. Thus, we may study the molec- 
ular structures of prion proteins to obtain some insights of prion diseases. Human 
prion diseases include the Creutzfeldt-Jakob disease, variant Creutzfeldt-Jakob dis- 
eases, Gerstmann-Straussler-Scheinker syndrome, Fatal Familial Insomnia, and Kuru. 
The NMR solution structure of the human prion protein (lQLX.pdb) was released into 
the PDB bank (www.pdb.org) in 1999 and last modified in 2009 |Zahn et al."2000] . 
Mice are popular experimental laboratory animals and the NMR structure of mouse 
prion protein (lAG2.pdb) was released into the PDB bank in 1997 and last modified in 
2009 [Ri ek et al. 1996] , Rabbits, dogs and horses were reported to be resistant to prion 
diseases [Vorberg et al. 2003} |Polymenidoua et al. 2008[ IPerez et al. 2010] and by the 



end of 2010 their NMR structures (2FJ3.pdb, lXYK.pdb, and 2KU4.pdb respectively) 
were completed to release into the PDB bank [Wen et al. 2010a] |Lysek et al. 2 005, 
IPerez et al. 2010] , The X-ray structure of rabbit prion protein (3079. pdb) was re- 
leased into PDB bank in 2010 too [K han et al. 2010] . At this moment it is very 
worth studying these molecular structures of horse, dog, rabbit, human and mouse 
prion proteins to reveal some secrets of prion diseases. The author found that dog 
and horse prion proteins have stable molecular dynamical structures whether under 
neutral or low pH environments | Zhang and Liu 2011 Zhang 2011a], but rabbit prion 



protein has stable molecular dynamical structures only under neutral pH environment 



Zhang and Liu 2011 Zhang 2011a Zhang 2010 Zhang 2011b . Under low pH envi 



ronment (at 450 K and 350 K), among all the prion proteins above-mentioned only for 
the rabbit prion protein, its stable a-helical molecular structures collapse into /3-sheet 
structures |Zhang and Liu 2011 [[Zhang 2011a] |Zhang 20101 |Zhang 20TTb"l . The conver- 



sion of disease PrP Sc from normal PrP is just involving 'conformational change' from 
predominantly a-helical protein to one rich in /3-sheet structure. This article specially 
focuses on the rabbit prion protein to obtain some insights into the immunity of rabbits 
to prion diseases. 

For the rabbit prion protein, we have its NMR and X-ray structures (2FJ3.pdb and 
3079. pdb respectively). Early in 2004, Epa |Zhang et al. 2006] made molecular model- 
ing of a homology structure (denoted 6EPA.pdb) for rabbit prion protein RaPrP (120- 
229), which was constructed using the NMR structure of the human prion protein 
(lQLX.pdb) as the template. Besides all these structures of rabbit prion protein, the 
knowledge on their conformational evolution/dynamics is considered essential to under- 
stand rabbit prion and the molecular modeling (MM) and molecular dynamics (MD) 
approach takes advantage beyond the experimental limit. In Section 2, this article first 



briefly reviews the main MD results of the homology MM structure |Zhang et al. 2 006 
at 500 K and of the wild-type NMR structure at 450 K |Zhang 201 lb| , compared with 
the MD of the wild-type NMR structures of human (HuPrP c (125-228), lQLX.pdb) 
and mouse (MoPrP c (124-226), lAG2.pdb) prion proteins. Section 3 will present the 
wild-type NMR rabbit, human and mouse MD comparisons at 350 K. Because the X- 
ray structure 3079. pdb was produced differently from the NMR structure 2FJ3.pdb 
and the Homology structure 6EPA.pdb, we will not do their MD comparisons; however, 
in Section 4, we will give detailed sequence and structure alignment analysis of all these 
three rabbit prion structures compared with human and mouse wild-type structures. 
Section 5 gives some concluding remarks on rabbit prion protein and prion diseases. 



MD Reviews On the 500 K Homology Rabbit Prion 
Protein and the 450 K NMR Rabbit Prion Protein 



2.1 The Homology Rabbit Prion Protein at 500 K [Zhang et aI7~2 006| 

Zhang et al. |Zhang et al. 2006| studied the MD of RaPrP c (120-229) homology struc- 
ture (6EPA.pdb). The MD simulations used Amber 8 [Case et al. 2004J PMEMD pro- 
gram, with explicit water at different temperatures and pH values. The simulation 
conditions are listed in Table [TJ The RMSD (root mean square deviation) and radius 



Table 1: Simulation conditions for the homology model. 



pH value 


Specie 


Truncated 


octahedral box (angstroms) Total atoms 


Ions added 


TIP3P Waters added 


Neutral pH 


HuPrP 






75.868 


19484 


3Na+ 


5929 




MoPrP 






67.447 


13422 


2Na+ 


3918 




RaPrP 






78.130 


21469 


2Na+ 


6572 


Low pH 


HuPrP 






74.834 


18530 


16C1- 


5599 




MoPrP 






67.335 


13208 


14U1- 


3836 




RaPrP 






80.896 


23847 


1401- 


7354 



of gyration results are shown in Fig. [TJ 

<Fi g . m> 

We may see that in Fig. [TJ rabbit prion protein has more stable structural dynamical 
behavior compared to the human and mouse prion proteins at 500 K under neutral pH 
environment. This is also shown in Fig. [2] of snapshots for human, mouse and rabbit 
prion proteins at 5ns, 10ns, 15ns, 20ns, 25ns, and 30ns respectively. 



<Fig. [2]> 



Fig. [2] shows that the helices of HuPrP and MoPrP were unfolded but RaPrP still 
keeps the helical structures at 500 K under neutral pH environment. Under low pH 
environment at 500 K, these helical structures of RaPrP were unfolded. One of the 
reasons of the rabbit prion protein unfolding is due to the remove of the salt bridges 
such as N177-R163 (Fig. EJ). 



<Fi g . G2> 

We may see in Fig. [3] the salt bridge / hydrogen bond between Arginine 163 and 
Aspartic acid 177 is conserved through a large part of the simulations and contributes 
to the protein stability of rabbit prion protein structure. Simulations at low pH value, 
where this salt bridge is absent, show RMSD and radius of gyration values for the rabbit 
prion protein to be of the same magnitude as the human and mouse prion proteins. 



2.2 The NMR Rabbit Prion Protein at 450 K [Zhang 2011b 



Zhang |Zhang 20Jl~b"1 did the MD studies on the NMR rabbit prion protein RaPrP c (124- 



228) (2FJ3.pdb) at 450 K under both neutral and low pH environments for the sim- 
ulations of 20 ns. Zhang |Zhang 201 lb| found that "the secondary structures under 
low pH environment at 450 K have great differences between rabbit prion protein and 
human and mouse prion proteins: the a-helices of rabbit prion protein were completely 
unfolded and began to turn into /3-sheets but those of human and mouse prion proteins 
were not changed very much. These results indicate the C-terminal region of RaPrP c 
has lower thermostability than that of HuPrP and MoPrP . Under the low pH envi- 
ronment, the salt bridges such as D177-R163, D201-R155 were removed (thus the free 
energies of the salt bridges changed the thermostability) so that the structure nearby 
the central helices 13 was changed for rabbit prion protein" |Zhang 2011b . The author 



continued his MD simulations for another 10 ns. The secondary structures for the MD 
simulations of 30 ns (Fig. [J]) shows the same conclusion as that of |Zhang 201 lb| . 

<Fig. 1> 

We may say that the salt bridges such as D177-R163, N201-R155 contribute to the 
structural stability of wild- type rabbit prion protein (Fig. U]). At 450 K, whether in 
neutral or in low pH environments, the a-helical secondary structures of dog prion 
protein have not changed for the long 30 ns' simulations [Zh ang and Liu 201 1| ; at 
350 K, horse prion protein has the same molecular structural dynamics |Zhang 2011a 
during the 30 ns' long simulations. 

3 350 K 

350 K might be a practical temperature for some experimental laboratory works. Zhang 
|Zhang 201 lc| did MD simulations for wild-type rabbit, dog and horse prion NMR struc- 
tures at 350 K. The findings of 350 K are: "dog and horse prion proteins have stable 
molecular structures whether under neutral or low pH environments. Rabbit prion pro- 
tein has been found having stable molecular structures under neutral pH environment, 
but without structural stability under low pH environment. Under low pH environ- 
ment, the salt bridges such as D177-R163 were broken and caused the collapse of the 
stable a-helical molecular structures" . Here the MD simulations are done for wild-type 
human and mouse prion proteins in the use of the same Materials and Methods as in 
Zhang 2011c . 



<Fig. E> 

Seeing Fig. [5j we know that at 350 K human and mouse prion proteins have stable 
molecular structures whether under neutral or low pH environments. 

Clearly the following salt bridges play an important role to the NMR structural sta- 
bility of rabbit prion protein: (1) GLU210-ARG207-GLU206-LYS203 (99.78%, 88.85%, 
82.74%, H3-H3), GLU210-HIS176 (74.31%, H3-H2), GLU206-HIS176 (57.10%, H3-H2), 
ARG207-HIS176 (0.52%, H3-H2), ASP177-ARG163 (19.54%, H2-S2); (2) ARG150- 
ASP146-ARG147-ASP143 (91.38%, 100%, 86.43%, Hl-Hl), HIS139-ARG150 (50.96%), 
HIS139-ASP146 (92.62%); (3) ASP201-ARG155 (10.07%, H3-H1), ASP201-ARG150 
(2.61%, H3-H1), ASP201-ARG147 (0.01%, H3-H1), ASP201-HIS186 (0.50%, H3-H2); 
and(4) ARG155-ASP201 (10.07%, H1-H3), TYR156-HIS186 (H1-H2, 71.69%), ARG155- 
GLU151 (20.70%, Hl-Hl), ARG155-GLU195 (0.06%), where HI, H2, H3 denote the 
a-helix 1, 2, 3 respectively, SI, S2 denote the /3-strand 1 and 2 respectively, and '%' 
denotes the percentage during the whole simulation of 30 ns. Compared with human, 
mouse, dog and horse NMR prion proteins, rabbit NMR prion protein has some special 
salt bridges which contribute to its structural stability at 350 K during the simulation 
of 30 ns (Fig. [6]) (human, mouse, dog and horse NMR prion proteins have not these 
salt bridges). 



<Fig. M> 



4 Alignment Analyses 



We make the sequence alignment of PrP from horse, dog, rabbit, human and mouse 
protein (Fig. [7|). 

<Fig. m> 

In Fig. "*" means that the residues in that column are identical in all sequences in 
the alignment, ":" means that conserved substitutions have been observed, "." means 
that semi-conserved substitutions are observed, the RED color takes place at small 
(small-|- hydrophobic (incl.aromatic-Y)) residues, the BLUE color takes place at acidic 
residues, the MAGENTA color takes place at Basic-H residues, GREEN color takes 
place at Hydroxyl+sulfhydryl+amine+G residues, and Grey color takes place at un- 
usual amino/imino acids etc.. For the structural domain, in Fig. [7] we can see some 
special residues listed in Table [2] for horse, dog, human and mouse prion proteins, which 
might contribute to characters of each structure respectively. Rabbits differ from horses, 
dogs, humans and mice at: S173 (N174 for horse, T174 for dogs, N174 for humans and 
mice), Q219 (K220 for horses and humans, R220 for dogs and mice), A224 (F225 for 
horses, Y225 for dogs, humans and mice), L232 (1233 for dogs, V233 for horse, humans 
and mice), and G228 (others are S229). For rabbits, at positions 89 and 97 the residues 
are special from all others (G89 (others are N90), S97 (others are N98)). These special 
residues are illuminated in Fig. [8j Some recent researches are focusing on the loop be- 



Table 2: Alignment analysis of special residues for HoPrP, DoPrP, HuPrP, and MoPrP. 

Horse S167 (others arc D), Y222 (others arc S), Q226 (others arc Y), V241 (others arc I), P245 (others arc S) 

Dog L129 (others arc M), S165 (others arc P), N170 (others arc S), S173 (others arc N), V244 (others arc I) 

Human 1138 (immunities arc L). S143 (others arc N), H155 (others arc Y), M166 (others arc V), 1183 (immunities arc V), 

E219 (others arc Q), S230 (immunities arc A) 

Mouse 1183 (immunities arc V), V215 (others arc I), D217 (others arc Q). S230 (immunities arc A) 



tween j32 and a2, i.e. PrP(164-171) |Apostol et al. 20114 IFernandez-Funez et al. 20TT] 



IKhan et al. "2()To| IPerez et al. "MlOJ IWen et al. 2"0l0a1 IWen et al. 2"oTOb] : we may see 
in Fig. [8] that the immune animals horses, dogs and rabbits have some residues in this 
loop different from humans and mice. 

<Fig. m> 

Lastly, we illuminate the figure (Fig. [9]) of rabbit prion protein, including the 
homology, NMR and X-ray structures (6EPA.pdb, 2FJ3.pdb, and 3079. pdb respec- 
tively). We superpose the homology structure onto the NMR structure and find the 
RMSD value is 3.2031669 angstroms. Similarly, we superpose the X-ray structure onto 
the NMR structure and we get their RMSD value is 2.7918559 angstroms. This im- 
plies the homology structure 6EPA.pdb made in 2004 by Epa |Zhang et al. 2006] is as 
effective as the X-ray structure 3079. pdb released recently on date 2010-11-24 (last 
modified on 2011-02-02). 



<Fig. [9]> 



5 Conclusion 



To really reveal the secrets of prion diseases is very hard. Prion proteins have two 
regions: unstructured region and structured region. Rabbits, horses, and dogs were re- 
ported having immunity to prion diseases. Fortunately, by the end of 2010 all the NMR 
molecular structures of rabbit, horse, and dog prion proteins had been released into 
PDB bank already; for rabbit prion protein, its X-ray structure was also released into 
PDB bank in the end of 2010. Prion diseases are 'structural conformational' diseases. 
This paper timely presents a clue to reveal some secrets in the view of the dynamics 
of prion molecular structures. MD experiences of the author nearly in the passing 10 
years show to us a common conclusion: under low pH environment at many levels of 
temperatures with different starting MD velocities, the rabbit prion protein always un- 
folds a-helical structures into /3-sheet structures. Prion diseases are just caused by the 
conversion from predominant a-helices of PrP into rich /3-sheets of PrP . Hence, we 
should furthermore study rabbits, horses and dogs, compared with humans and mice 
in order to reveal some secrets of prion diseases; for us, it is a long shot but certainly 
worth pursuing. 
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Figure 1: RMSD and Radius Of Gyration graphs at 500 K neutral pH value. 
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Figure 2: Snapshots of HuPrP, MoPrP and RaPrP (from up to down) at 500 K neutral 
pH value, at 5ns, 10ns, 15ns, 20ns, 25ns, and 30ns respectively (from left to right). 
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Figure 3: Distance between D177.0D1 and R163.NE for RaPrP, between D178.0D1 
and R164.NE for HuPrP and MoPrP, at 500 K neutral pH value. 
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Figure 4: Secondary structures of rabbit, human and mouse prion proteins (from up 
to down) at 450 K under neutral and low (from left to right) pH environments (red: 
a-helix, pink: ir- helix, yellow: 3io-helix, green: /3-bridge, blue: /3-sheet, purple: Turn, 
Black: Bend; x-axis: time (0-30 ns), y-axis: residue numbers). 
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Figure 5: Secondary structures of rabbit, human and mouse prion proteins (from up 
to down) at 350 K under neutral to low pH environments (from left to right) (X-axis: 
ns - 30 ns (from left to right), Y-axis: residue numbers 124 - 228 / 125 - 228 / 124 - 
226 (from up to down)). 
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Figure 6: Some special salt bridges ARG207-HIS176, TYR156-HIS186, HIS139- 
ARG150, ASP201-ARG147, ASP201-ARG150, ASP201-HIS186, ARG155-GLU151 of 
wild-type NMR rabbit prion protein at 350 K. 



15 



CLUSTAX 2 . 1 noltiple frequence alignment 

&0123456Ta&Gl2345H$?a3<H234567flSO L234567B33U2345679901234561 
Horse PrP MVKSRVGGWrLVI,FVATW3I^GLCtUtRPKPCG-i™TGG3RTP&QCSPGeWRYPPQGeGGW 55 

DogPrP MVKSKIGSWILV[»FVAMWSDVGLCKKRPKPGC;GWWTGGSR,i!Pi3QGSPC?GNRYPPQGGGGH 60 

RiibbiLPi-P — MAHLGYW^IAWATWSWGLCKKRPKFGGGWlTrGGSRYPGQSSPGGWRYPPQGGG-W 51 

HUtnanPtP MANIiGCWMLVl.FVftTWaDLiGLCKFRPKPGG WNTGGSRYPGQGSPGGNRYPPQGGGGW 57 

MOUSSPrP - -MANLGYWLLAI.ETVTMWTDVGtCElKRPKPGG-WHTGGSRYPGQGSPGGMRYPPgGG-TH 56 

; . . * *.# *#*. * -y . *** ******** ***********_***********# * 

HoraePrP G3PRGGGWGQPHGGGWG£PHGGGtoGQPHGGGGWGGGG- SHGOWNKPSKPrcrNHKJIVAGAA 118 

DogPrP GQPHGGCTfGQPHGGGWG^HGGGWG^PHGGGGWGiMG-THSQWWKPSKPKraMKllVAGJiA 119 

Rabb it P r P GQPHGGGWGQPHGGGWGQPHGGGWGQPHGGG - WGQGG THNQWGK PSK PKTSMKITVAGAA 1 IS 

Human PrP GQPHGGGWGQPHGGGWGQPHGGGHGQPHGGG HGQGGGTH3QWNKPSKPKTNMKKMAGAA 116 

Mouse PrP GQPHGGGHGQPHGGSWGQFHGGSWGQPHGGG WG^GGTHNQWKKPSKPKTNLKHVAGAA 115 

+**+*********+ ******* _* + * ***** ***** -* i ±* i ***i+*+ i -+*-**** 

119 130 140 150 160 170 

901234 56789012 34567B9Q123 456789Q12 34 56 78901234567B901 2 3456 

HorsaPrP AAGAWGGLGGYl€LG5M43FJLIllFGNDYEDRYYRET511OTRYFNQVTYRFV5EY5flQNNFV 173 

DogPrP AAGA WGGLGGY I.LG3 AW3 RFL I ] iFGN DCEDRY YKENMYRY FNQV Y YH 3 V DQ YNNQ3 G' FV 179 

RabbitPrP AAGAWGGLGGYMLGSflMSSPLIliFGblDYEDRYYHENMYRYPNQVYIfUPVDQYSNgiMSFV 175 \ 

HtimanPrP AACAWGGLGGYMI.G3AM3RP1 IHFGSDYEDiWfYRENMHRYPMQVY YRPMDEY3NQWNFV 17 6 

MousePrP AAGAVVGGLGGYMLSSAHSHPMIHFGHDMEDRYYRENMYRYPHQVYYRPVDQYSNQNNPV 175 

**A***-***.-* H*. • **i***#i .* * + * * *#*•##**»■*.*■****** + * - ■* *+ *# 



HorsePrP 
DogPrP 
RabbitPrP 
Human PrP 
MousePrP 



LflO 190 20D 210 220 231 

78 9012345678 901234 56739012 34 56783012 345678901234 567 3 90133-1 



HDCVN I TVK0HT VTTTT KGENFTETDVK I MERWEQMCI TQ YQEE YEAFqQ - HG AS WL 236 
HDCVN I TVKQflT VT'FT- KGENFTETDIKMMERWEQMCI TQY0RESER YYQ --RGASVTL 236 
HDCVN 1 1' VKftliT V1TT TXGENFTIJTD IK LMEKWBQMC 1 TQ YQQESQAA YQ - - RAAG V1.L 233 
I1DCVN I T IKQHT VTTTTKGENFTETD VKMMERWEQNCI TQYERE3QA YYQ — RGS SMVL 234 
HDCVN 11' IKQH1 1 VTTl'TKGBNFTETDVKMHERWE^MCVTQYQKEBQAY Y DGRRSS ST VL 235 

tH**«l.M+*4*H a********.*.*********.***..* ■« ■ * . -4 



Horse PEP 

DogPrP 

RabbitPrP 

HumanPrP 

MousePrP 



'■■ ','fi'Sitl.v i'1'i<>7r;"in> I 

F3SPPWLLIFFLIFLIVG 255 

FSSPPVJI.LVSFLTFLIVG 255 

FSSFFVILLISFLIFLIVG 252 

F3SPPVIl.LT3FLIFI.IVG 253 

FSSFFVILLISFLIFLIVG 254 
******-**- ******** 



Figure 7: Horse, dog, rabbit, human and mouse prion protein sequence alignment. 
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Figure 8: Special residues owned only by HoPrP, DoPrP, RaPrP, HuPrP, and MoPrP 
respectively. 
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Blue: NMR Structure 
Green: Homology Structure 
Red: X-ray Structure 




Rabbit PrP 



Figure 9: Rabbit prion protein NMR, homology and X-ray structures (2FJ3.pdb, 
6EPA.pdb, and 3079.pdb). 
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